Roughening transition driven by a binary spreading process 
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We introduce a solid-on-solid growth process which evolves by random deposition of dimers, 
surface diffusion, and evaporation of monomers from the edges of plateaus. It is shown that the 
model exhibits a robust transition from a smooth to a rough phase. The roughening transition is 
driven by an absorbing phase transition at the bottom layer of the interface, which displays the 
same type of critical behavior as the pair contact process with diffusion 2A — > 3A, 2A — > 0. 
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I. INTRODUCTION 

The fabrication of thin films by deposition and evapo- 
ration of particles plays an important role in technologi- 
cal applications such as molecular beam epitaxy. Theo- 
retical studies of growing films are usually based on sim- 
plified models of interfaces evolving by certain stochas- 
tic dynamic rules. Even simple solid-on-solid models 
may display very interesting features such as kinetic 
roughening, scaling, self-similarity, and fractal proper- 
ties |, |,||, | . 

Particularly fascinating are models with a roughen- 
ing transition from a smooth phase, where the width of 
the interface is stationary and finite, to a rough phase, 
where the width grows continuously. Roughening tran- 
sitions may rely on different physical mechanisms. For 
example, Kardar-Parisi-Zhang (KPZ) growth processes 
are known to exhibit a roughening transition in d > 2 
dimensions, which is caused by a competition of surface 
tension and the nonlinear term in the equation of motion. 
Another class of roughening transitions, on which we will 
focus in the present work, occurs in certain deposition- 
evaporation processes with the special feature that atoms 
cannot desorb from the middle of plateaus. In these mod- 
els a completed layer, once formed, is stable and cannot 
evaporate again. Thus, if the growth rate is sufficiently 
small, the interface is smooth and fluctuates close to the 
lowest completed layer (the so-called bottom layer of the 
interface) , while for large deposition rate the growing is- 
lands coalesce so that the interface propagates and rough- 
ens continuously. 

A very simple model, which exhibits such a roughening 
transition, has been introduced some time ago by Alon et 
al. U . The model is defined as a restricted solid-on-solid 
growth process with the special property that atoms can 
only evaporate only at the edges of plateaus. It turned 
out that the roughening transition in this model is driven 
by a directed percolation (DP) process at the bottom 
layer. DP is the standard universality class of phase tran- 
sitions from a fluctuating phase into an absorbing state 
and occurs in unary reaction diffusion processes of the 
form A -> 2 A, A -> |. If the decay A -> dom- 
inates the system eventually enters an absorbing state 
without particles from where it cannot escape, while in 



the active phase the reaction A — > 2A is strong enough 
to sustain a stationary density of particles. As shown in 
Ref. [H the sites at the bottom layer of the interface can 
be interpreted as particles of a DP process. Using this 
relation the pinned phase, where the interface is smooth, 
corresponds to the active phase of DP. However, if the 
growth rate is increased above a certain critical thresh- 
old the DP process eventually enters the absorbing state, 
meaning that the bottom layer propagates by one step. 

The relation to DP allows one to predict some of the 
critical properties at the roughening transition. Since the 
density of sites at the bottom layer no should scale in the 
same way as the particle density in DP, it is expected to 
decay at criticality as 

n ~ t~ 5 , (1) 

where 5 = 0/v± is one of the universal critical exponent 
of DP. Extending these arguments it has been shown that 
even the critical behavior at the first few layers above 
the bottom layer can be described in terms of a unidirec- 
tionally coupled hierarchy of DP processes M . Roughly 
speaking each layer is associated with a separate DP pro- 
cess which is coupled to the dynamics at the layer below, 
while effective couplings in opposite direction turn out to 
be irrelevant in the renormalization group sense. 

More recently it became apparent that other universal- 
ity classes of phase transitions into absorbing states can 
also be used to generate a roughening transition in an 
appropriately defined growth model. For example, the 
so-called parity-conserving (PC) class M, which is most 
prominently represented by branching-annihilating ran- 
dom walks with two offspring A — > 3A, 2A — > 0, has been 
shown to generate a roughening transition in a simple 
model for dimer adsorption and desorption [lC| ]. There- 
fore, the question arises whether it is possible to relate 
any universality class of absorbing phase transitions to a 
corresponding class of roughening transitions. 

As a complete classification scheme of absorbing phase 
transitions is not yet available, it is an important task 
to search for other yet unknown types of transitions. In 
this context one of the most promising and controver- 
sially debated candidates is the pair contact process with 
diffusion (PCPD) 

2A ->3A, 2A -> . (2) 



(a) Deposition 



(b) Evaporation 



(c) Hopping 



FIG. 1: List of all possible rearrangements of the interface at the sites i — 1, i, i + 1, i + 2 according to the dynamic rules ([J)-(|7|): 
(a) Deposition of a dimer at the sites i and i + 1. (b) Evaporation of a monomer from the right edge of a terrace at site i. 
(c) Hopping of a monomer between two adjacent islands, moving the hole between the islands in opposite direction. In all 
cases the spatially reflected rules (not shown here) have to be included as well. 



In contrast to DP and PC transitions, where individual 
particles can generate offspring, the PCPD is a binary 
spreading process, i.e., two particles have to meet in or- 
der to create a new one. Although there is a still ongoing 
debate to what extent the PCPD represents a genuine in- 
dependent universality class, there seems to be a general 
consensus that it displays a novel type of phase transi- 
tion that has not been seen in other models ITT[ |l2| , |l5| , 
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numerically accessible temporal regime the transition is 
characterized by a fairly reproducible set of effective crit- 
ical exponents. Very recently Kockelkoren and Chate ||f| 
carried out extensive simulations, reporting the estimates 



= 0.37(2), 5 = — = 0.200(5), z = ^- = 1.70(5) . 

(3) 

The purpose of this paper is to introduce a solid-on-solid 
growth process which exhibits a roughening transition 
driven by a binary spreading process at the bottom layer. 
Focusing on the one-dimensional case, where fluctuation 
effects are most pronounced, we demonstrate that the 
pair contact process with diffusion can indeed be used 
to generate a roughening transition in an appropriately 
defined model for interface growth. 



II. DEFINITION OF THE MODEL 

The dynamic rules of the model involve three differ- 
ent physical processes, namely, deposition of dimers, sur- 
face diffusion, and evaporation of diffusing monomers. 
The dimers are deposited horizontally on pairs of sites at 
equal height, leading to the formation of islands. These 
islands arc stable in thc interior but slightly unstable 
at the edges, where monomers are released at a certain 
rate. These monomers then diffuse on the surface until 
they either attach to another island or evaporate back 
into the gas phase. In what follows we consider the limit 



of a very high rate for evaporation. In this limit a re- 
leased monomer is most likely to evaporate unless it im- 
mediately attaches to an adjacent island at the next site, 
effectively moving the hole between the two islands in 
opposite direction. 

In more technical terms, the model is defined on a 
one-dimensional lattice with periodic boundary condi- 
tions. Each site i is associated with an integer height 
variable hi . The model evolves by random-sequential up- 
dates, i.e., a site i is randomly selected and one of the 
following moves is carried out (see Fig. |l|): 

(a) With probability q a dimer is deposited horizon- 
tally provided that the two supporting sites have 
the same height: 



If hi — hi + i then 
hi — *■ hi + 1; fti+i 



H+l 



1 



(4) 



(b) With probability monomer is released at the 

edge of a terrace. If there is no adjacent terrace, 
the released monomer evaporates immediately: 



If hi > hi+i > hi + 2 then 
hi—>hi — l 

If hi + \ > hi > hi^i then 
hi+i — » hi + \ — 1 



(5) 



(c) Otherwise the released monomer attaches to the 
opposite edge of the adjacent terrace: 



If hi > hi+i < hi + 2 then 

hi -> hi - 1; hi+i -> hi+i + 1 

If hi+i > hi < hi— % then 

hi+i — > hi+i — 1; hi — > hi + 1 



(6) 



In addition, an attempted move is abandoned if the re- 
sulting configuration would violate the restricted solid- 



q=0.75 

FIG. 2: Typical interface configurations for various values of 
the growth rate q. 

on-solid (RSOS) condition 

\hi -h i±1 \ < 1. (7) 

The RSOS condition introduces an effective surface ten- 
sion and imposes an infinitely large Ehrlich-Schwoebel 
barrier, i.e., particles cannot diffuse across the edges of 
terraces. Each attempted update corresponds to a time 
increment of At = 1/L. Note that the dynamic rules 
are translationally invariant in horizontal as well as in 
vertical direction. 



III. PHENOMENOLOGICAL PROPERTIES 
AND RELATION TO THE PCPD 

The growth process defined above has the following 
phenomenological properties (see Fig. |^). For small val- 
ues of q occasionally deposited dimers will quickly desorb 
from the surface by evaporation at the edges so that the 
interface is smooth and fluctuates close to the bottom 
layer. As q is increased the islands grow and coalesce 
until their average size diverges at a certain critical point 
q c . Above the critical point more and more islands are 
deposited on top of each other so that the interface rough- 
ens continuously. However, the interface remains pinned 
at solitary sites at the bottom layer, leading to the forma- 
tion of pyramidial structures. These pyramids fluctuate 
so that the pinning sites slowly diffuse until they merge 
or annihilate. In contrast to the dimer model introduced 
in Ref. jf(| there is no evidence for a facetting transition 
at q = 1/2. 

In order to explain how the growth model is related to 
the PCPD, let us consider the dynamic processes at the 
bottom layer (the spontaneously selected lowest height 
level of the interface). Interpreting sites at the bottom 
layer as particles A of a spreading process, the dynamic 
rules listed in Fig. [I] can be associated with certain reac- 
tions of the particles. For example, the deposition of a 
dimer corresponds to pairwise annihilation 2 A — > 0, while 
evaporation of a monomer can be viewed as the creation 
of a particle A. However, in the present model atoms 
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FIG. 3: Density of sites at the lowest height (bottom layer) as 
a function of time for q = 0.15900, 0.15910, 0.15915, 0.15920, 
0.15925, 0.15930, and 0.15940 averaged over up to 1000 inde- 
pendent runs. The inset shows the same data multiplied with 
t s , where <5 = 0.21(2) is the estimated exponent. 

can only evaporate at the edge of a terrace followed by 
two bottom layer sites, hence rule (b) in Fig. |I] has to be 
interpreted as a fission process 2A — > 3A Otherwise, if 
there is only one such bottom layer site next to the edge, 
rule (c) applies which corresponds to a random walk of 
a single particle A. Thus we may interpret the dynamics 
at the bottom layer as a reaction-diffusion process 

(a) AA -> 00 

(b) AA% I %AA -> AAA (8) 

(c) 0A «-» A% 

which resembles the dynamic rules of the pair con- 
tact process with diffusion in a fermionic realization. 
Clearly, the relation between the growth model and the 
PCPD is not rigorous, especially because of deposition- 
evaporation processes on top of islands. However, as in 
the case of DP- and PC-related growth models |j, [To| |, 
this correspondence is expected to be valid in the renor- 
malization group sense and therefore determining the 
asymptotic critical behavior. 



IV. CRITICAL PROPERTIES AT THE 
ROUGHENING TRANSITION 

A. Scaling at the critical point 

In order to confirm that the roughening transition in 
the present model is driven by a binary spreading pro- 
cess, we show that the density of sites at the bottom 
layer no displays the same type of decay as the density of 
particles in the PCPD. To this end we perform numerical 
simulations starting with an initially flat interface. Mea- 
suring no for various values of q, the critical threshold q c , 
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FIG. 4: Mean height H(t) and squared width W 2 (t) as func- 
tions of time in a semilogarithmic plot for the same parame- 
ters as in Fig. []. 



where the density decays approximately as a power-law, 
is estimated by 



q c = 0.15920(5). 
The corresponding exponent 

5= — = 0.21(2) 



(9) 



(10) 



is indeed in good agreement with the expected PCPD 
exponent (||). However, the plotted lines are still slightly 
curved, a problem which also occurs in many other binary 
spreading processes. 

Fig|| shows the the average height 



#(t) = 7$>(t) 



(11) 



and the squared width 



W 2 {t) 



ij2Ht)-H(t)) 2 (12) 



as functions of time. In previously investigated DP- and 
PC-driven roughening transitions these quantities were 
found to increase logarithmically with time at the critical 
point. In the present case the average height shows an al- 
most logarithmic increase H (t) ~ 0.19+0.13 lni, whereas 
the squared width W 2 (t) does not display a convincing 
logarithmic law, at least within the numerically accessi- 
ble range. 
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FIG. 5: Finite size scaling. Data collapses according to the 
scaling forms (|l|), (|l|) and (0) for system sizes L=32, 64, 
128, . . ., 1024 averaged over 15000 independent runs. Data 
points for t < 50 are eliminated in order to suppress deviations 
due to initial transients. Moreover, the curves are cut for 
t > 4 x L 1 ' 76 , from where on the interface may become pinned 
at a single bottom layer site, leading to additional deviations 
for large t. 



scaling theory the bottom layer density no is expected to 
obey the finite-size scaling form 



n (L,t) = t- 5 f(t/L z 



(13) 



where f(z) is a universal scaling function. This scaling 
form can be used to generate a data collapse by potting 
n (L, t)t s versus t/L z , as shown in the left panel of Fig. ||. 
Using the previous estimate (|l0|) and tuning z we obtain 
an excellent data collapse for 



z = 1.75(2) . 



(14) 



which is in agreement with the value in Eq. (||) for the 
PCPD. As in DP- and PC-driven roughening transi- 
tions, we expect this exponent to be valid not only at 
the bottom layer but for the growth process as a whole. 
Therefore, the dynamic exponent z should also determine 
the finite-size scaling properties of the squared width 
W 2 (L,t) and the average height H(L,t). In the case of 
DP-driven roughening transitions these quantities were 
found to obey the scaling forms p7| 



exp(W 2 {L,t)) 
exp(H(L,t)) 



L x Fi(t/L z 
L»F 2 (t/L z 



(15) 
(16) 



so that it is natural to assume that the same forms hold 
in the present case. In fact, plotting W 2 (L, t) — A In L and 
H(L,t) — /ilnL versus \n(t/L z ) both quantities show a 
reasonable data collapse for 



B. Finite-size scaling 

In order to estimate the dynamic exponent z we per- 
formed finite-size simulations. According to the standard 



A = 0.08(2), n = 0.21(3) 



(17) 



as shown in the right panel of Fig. ||. However, the col- 
lapse is less precise as the one for no because of severe 
corrections to scaling, as mentioned before in Fig. m 
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FIG. 6: Off-critical simulations. Data collapses according 
to the scaling forms @-© for e = ±0.0004, . . . , ±0.0128. 
Data points for t < 25 are eliminated in order to avoid devi- 
ations due to initial transients. 

C. Off-critical simulations 

So far we determined two ratios of the three standard 
exponents (3, i/±, vn. In order to determine the third ex- 
ponent one has to perform off-critical simulations. Here 
we expect the scaling forms Q 

n (e,t) = t- s g(t\ep) (18) 
exp(W/ 2 (e,i)) = lel-^GxMep) (19) 
exp(iJ(e,i)) = \e\-^G 2 (t\eP), (20) 

where e = q — q c denotes the distance from criticality. 
These scaling forms should hold below and above criti- 
cality, at each case with different scaling functions. The 
corresponding data collapses are shown in Fig. |[ It turns 
out that the best collapses for the bottom layer density 
no are obtained for 

2/|| = 1.85(5). (21) 

In the smooth phase (lower branches) the interface width 
and the average height show reasonable data collapses 
for \v ± = 0.09(2) and uv± = 0.235(15). With v ± m 1.1 
(taken from Eq. (^)) these estimates are compatible with 
Eq. (p7[). In the rough phase (upper branches) there is 
no convincing data collapse, indicating that the roughen- 
ing interface crosses over to a different type of dynamic 
critical behavior, probably related to the KPZ equation. 



V. CONCLUSIONS 

To summarize we have introduced a solid-on-solid 
growth process which exhibits a novel type of roughen- 
ing transition in 1+1 dimensions. The dynamic rules are 
defined in such a way that the sites at the bottom layer 
(the spontaneously selected lowest layer of the interface) 
evolve effectively as a pair contact process with diffusion 
2A — > 3^4, 2A — > 0. Measruring the bottom layer density 



no in finite-size and off-critical Monte Carlo simulations 
we have estimated the critical exponents 

(5 = 0.21(2), z= 1.75(2), ^= 1.85(5). (22) 

These exponents are in good agreement with recent esti- 
mates for the PCPD reported in |25), confirming that the 
roughening transition is in fact driven by a pair conact 
process at the bottom layer. By comparing these expo- 
nents we implicitly assume that the PCPD is character- 
ized by simple power-law scaling with well-defined critical 
exponents, a postulate which is still debated jlj| Q p5| . 
But even in terms of effective exponents the numerical 
results clearly confirm that the roughening transition in 
the proposed model is in fact driven by an underlying 
pair contact process with diffusion. 

In the active phase of the PCPD the interface is smooth 
an pinned to the bottom layer, while it roughens in the in- 
active phase. Since in the inactive phase of the PCPD the 
particle density decays algebraically as no ~ 1/%/t, the in- 
terface does not propagate uniformly, instead it remains 
pinned and becomes increasingly facetted. At criticality 
the interface width and the average height of the interface 
grow logarithmically with time or even slower. Adopting 
the scaling picture outlined in Ref. p7[ , the exponenti- 
ated quantities exp[PF 2 (i)] and exp[H(t)] are expected to 
obey ordinary power-law scaling with corresponding crit- 
ical exponents A ~ 0.08 and a ~ 0.21, respectively. How- 
ever, these estimates have to be taken with care as both 
quantities, especially the interface width, show strong de- 
viations from power-law scaling. 

Unlike the dimer model introduced in jlFj] , the present 
model does not exhibit a facetting transition. In the PC- 
driven case, where dimers are deposited and evaporated, 
the facetting transition is related to a biased diffusion 
of steps along inclined parts of the interface with slope 
1. Such a step moves upward (downward) if a dimer is 
deposited (evaporated) so that for q > 1/2 all terraces are 
swept to the top of the pyramids, stabilizing the facetted 
state. In the present model, where monomers evaporate 
instead of dimers, there is no such mechanism. 

As in the case of DP- and PC-driven growth processes, 
the critical behavior at the first few layers above the bot- 
tom layer may be described by a unidirectionally cou- 
pled sequence of diffusive pair contact processes. In fact, 
the densities ni(t) and n%{t) (not shown here) display 
the same qualitative behavior as the particle densities in 
such a coupled hierarchy. However, a systematic study 
of a coupled sequence of PCPD's should be postponed 
until the PCPD itself is fully understood. 

The results of the present work suggest that any ab- 
sorbing phase transition can be used generate a rough- 
ening transition in an appropriately defined growth pro- 
cess. For example, a generalization of the present model 
to a recently introduced triplet process Q is straight 
forward. Moreover, it would be interesting to study the 
influence of evaporation from the middle of plateaus com- 
bined with a hard-core wall at zero height, as it was done 
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in the case of DP-related growth processes in Refs. p9|| . 
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